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Problematic behaviours are severe welfare issues for one of the world’s most popular pets, the domestic 
dog. One of the most prevalent behavioural problem that causes distress to dogs is social fearfulness, 
meaning fear of conspecifics or unfamiliar people. To identify demographic and environmental factors 
associated with fear of dogs and strangers, logistic regression was utilised with a large dataset of 

6,000 pet dogs collected through an owner-filled behavioural survey. Social fearfulness was associated 
with several factors, including urban environment, poor socialisation during puppyhood, infrequent 
participation in training and other activities, small body size, female sex, and neutering. In addition, 

we identified several breed differences, suggesting a genetic contribution to social fearfulness. These 
findings highlight the role of inadequate socialisation, inactivity, and urban living environmental in 
fear-related behavioural problems in dogs. Improvements in the management and breeding practices of 
dogs could, therefore, enhance the welfare of man’s best friend. 


Fear is a major welfare problem in pet dogs’. As a transient feeling aroused by specific stimuli, fear is a normal, 
fundamental emotion conserved among species which aids an individual to survive from threatening situation”. 
Fearfulness, on the other hand, is a personality trait®. If fear is excessive, prolonged, or generalised in nature, 
fearfulness becomes a behavioural problem which can interfere with the normal performance of the dog, causing 
high levels of distress or anxiety and increasing the risk of diseases or even decrease lifespan’®. Moreover, prob- 
lematic behaviours might have a negative valence not only on the wellbeing of the dog but also on the wellbeing 
of its owner®°. Undesirable behaviours, such as excessive fearfulness, are the leading cause for relinquishment 
or even euthanasia of pet dogs worldwide!*. In the worst case, fearful dogs may resort to aggression and cause 
public health threats!>~”’, highlighting the gravity of canine fearfulness. 

Based on the stimulus eliciting fear, fearfulness can be divided into two separate categories in dogs: social 
and non-social fearfulness°. Social fearfulness encompasses fear of conspecifics or unfamiliar people, whereas 
non-social fearfulness includes fear of different stimuli, such as loud noises, novel situations, or heights and sur- 
faces. Fearfulness is one of the most frequent canine behavioural problems as the prevalence of fearfulness ranges 
from 26.2% even to 44%'*!?, and around 10-19% of dogs show fear of strangers or dogs'”!*°. 

Many behavioural traits are complex. Multiple genes with only small effects”!, a plethora of environmen- 
tal factors”, and the complex interplay between these genetic and non-genetic factors” all contribute to the 
development behaviour. Fearfulness in dogs is moderately heritable and associated with some candidate loci and 
genes**** in addition to environmental effects”””?-*’. However, more research is needed to reveal the genetic risk 
variants and environmental factors associated with social fearfulness in dogs. Enhanced understanding of these 
genetic and environmental factors interacting together, leading to the onset of a behavioural problem, would give 
us tools to better recognise, manage, and prevent these conditions. 

Asa part ofa larger population-based canine behavioural survey with over 13,700 participants, the aim of this 
study was to investigate the demographic and environmental factors associated with social fearfulness in Finnish 
pet dogs. By identifying the factors increasing the risk for social fear-related problems, we could improve the 
wellbeing of pet dogs. 
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Fear of dogs Fear of strangers 
Variable xX? p-value DF x? p-value DF 
Age 8.40 0.004* 1 0.24 0.626* 1 
Age\2 12.67 <0.001* 1 3:27 0.070* 1 
Sex 31.75 <0.001* 1 4.86 0.027* 1 
Sterilisation 48.65 <0.001 1 29.52 <0.001 
Sex*Sterilisation 12.29 0.003 1 
Breed 128.30 <0.001 22 131.57 <0.001 9 
Body size 115.96 <0.001 2 10.15 0.027 2 
Socialisation score 82.91 <0.001* 1 129.71 <0.001* 
Urban environment score 33.46 <0.001 1 14.52 0.001 
Activities/training 29.20 <0.001 2 15.39 0.003 2 
Daily exercise 13.36 0.015 3 
Family size 11.81 0.067 4 
Weaning age 8.93 0.097 3 


Table 1. Associations between the demographic and environmental variables with fear of dogs and fear of 
strangers in the logistic regression analyses. P-values are controlled for false discovery rate except for a priori 
contrasts. A priori effects are denoted with. *Significant effects are emboldened (p-value < 0.05). N = 5,973 
(fear of dogs) and N = 5,932 (fear of strangers). 


Results 

Study cohort and demographics. We studied the demographic and environmental factors associated 
with fear of dogs and fear of strangers in datasets of 5,973 and 5,932 dogs, respectively. In the ‘fear of dogs’ data, 
the numbers of non-fearful and fearful dogs were 4,806 and 1,167, respectively. The age of the dogs in this dataset 
varied from 2 months to 17 years (mean age 4.6 years). In the ‘fear of strangers’ data, the numbers of non-fearful 
and fearful dogs were 5,036 and 896, respectively. The age of the dogs in this dataset varied from 2 months to 17 
years (mean age 4.7 years). 51% of the dogs were females in both datasets. More detailed demographics and the 
lists of included breeds and the number of individuals per breed are presented in the Supplementary Table S1. 


Demographic and environmental factors associated with fear of dogs. Logistic regression analysis 
identified several demographic and environmental factors associated with fear of dogs, including socialisation 
score, breed, body size, age, urban environment score, activities/training, daily exercise, and the interaction of sex 
and sterilisation (Table 1, Fig. 1). 

In the logistic regression model, socialisation score, body size, and breed had the strongest associations with 
fear of dogs. Dogs with less socialisation during puppyhood were more likely to show fear of dogs (x* = 82.91, 
DF = 1, p < 0.001) (Table 1, Fig. 1). Small dogs were more likely fearful when compared to both large (OR 
= 3.29, p < 0.001) and medium sized (OR = 2.04, p < 0.001) dogs, and there was also a significant differ- 
ence between medium and large dogs (OR = 1.61, p < 0.001) (Supplementary Table S2, Fig. 1). We discovered 
behavioural differences between breeds. Chihuahua, Shetland Sheepdog, and Spanish Water Dog were the most 
fearful breeds whereas Pembroke Welsh Corgi, Cairn Terrier, and Wheaten Terrier were the least fearful breeds 
(Fig. 1). The largest pairwise differences were seen between Chihuahua and Pembroke Welsh Corgi (OR = 8.64, 
p < 0.001), Shetland Sheepdog and Pembroke Welsh Corgi (OR = 8.13, p < 0.001), and Spanish Water Dog 
and Pembroke Welsh Corgi (OR = 7.32, p < 0.001). Cairn Terrier also had a significantly lower likelihood of 
being fearful when compared to Chihuahua (OR = 0.13, p < 0.001), Shetland Sheepdog (OR = 0.14, p < 0.001), 
and Spanish Water Dog (OR = 0.15, p < 0.001). Significant pairwise breed differences are summarised in the 
Supplementary Table S3 and all pairwise breed differences are presented in the Supplementary Dataset. Based on 
previous research on breed differences of social fearfulness, we made an a priori hypothesis that Chihuahua, Jack 
Russell Terrier, Lagotto Romagnolo, and Shetland Sheepdog would be more fearful than German Shepherd Dog, 
Golden Retriever, Labrador Retriever, and Staff. Bull Terrier. Indeed, the first breed group was more likely afraid 
of dogs than the latter breed group (OR = 1.57, p = 0.024) (Supplementary Table S2). 

There was an association between the age of the dog and fear of dogs, indicating that dogs from two to eight 
years old had the highest probability of being fearful, but this likelihood decreased after eight years of age (linear 
effect: x* = 8.40, DF = 1, p = 0.004; quadratic effect: x? = 12.67, DF = 1, p < 0.001) (Table 1, Supplementary 
Fig. S1). There was a significant interaction between sex and sterilisation of the dog, as intact males were less likely 
to show fear of dogs when compared to intact females (OR = 0.62, p < 0.001), but no difference was observed 
between neutered males and females (OR = 1.10, p = 0.807). Moreover, intact dogs were less likely fearful than 
neutered dogs in both sexes (intact male vs. neutered male: OR = 0.45, p < 0.001; intact female vs. neutered 
female: OR = 0.76, p = 0.036) (Supplementary Table S2, Fig. 1). 

Dogs living in a more urban environment (x* = 33.46, DF = 1, p < 0.001) (Table 1, Fig. 1), participating less 
frequently in activities and training, and getting less daily exercise had higher probabilities of being afraid of dogs 
(Table 1, Supplementary Table S2, Supplementary Fig. $1). Dogs participating in activities only seldom or never 
were more likely fearful than dogs training sometimes (OR = 1.54, p < 0.001) or weekly (OR = 1.56, p < 0.001). 
Dogs getting less than one hour of exercise per day were more likely fearful than dogs exercising more than three 
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Figure 1. The effects of breed, socialisation, interaction between sex and sterilisation, body size, and living 
environment on fear of dogs in the logistic regression analysis. (a) Breeds differed in the likelihood of showing 
fear of dogs, with Chihuahua being the most and Pembroke Welsh Corgi the least fearful breed. (b) Dogs that 
had less socialisation experiences at the age of 7-16 weeks were more afraid of other dogs. (c) Intact males 
were less likely afraid of dogs than intact females, but no difference was observed between neutered males and 
females. In addition, intact individuals were less fearful in both sexes when compared to neutered individuals. 
(d) Small sized dogs were more likely fearful than medium and large sized dogs. There was also a difference 
between medium and large sized dogs. (e) Dogs living in a more urban environment had a higher likelihood of 
showing fear of dogs. Grey lines (b,e) and error bars (a,c,d) indicate 95% confidence limits. N = 5,973. 


hours daily (OR = 1.58, p = 0.022). Dogs exercising 1-2 hours per day were also more likely fearful when com- 
pared with dogs exercising more than three hours in a daily basis (OR = 1.50, p = 0.005). 


Demographic and environmental factors associated with fear of strangers. The best model 
explaining the difference between fearful and non-fearful dogs included several demographic and environmental 
factors, including socialisation score, breed, age, sex, sterilisation, body size, urban environment score, activities/ 
training, weaning age, and family size (Table 1, Fig. 2). 

Socialisation score and breed had the strongest associations with fear of strangers in the logistic regression 
analysis (Table 1, Fig. 2). Dogs with less socialisation experiences were more likely afraid of strangers (x? = 
129.71, DF = 1, p < 0.001) (Table 1, Fig. 2). We discovered behavioural differences between breeds. Spanish 
Water Dog, Shetland Sheepdog, and Chinese Crested Dog were the most fearful breeds whereas Wheaten Terrier, 
Finnish Lapponian Dog, and Labrador Retriever were the least fearful breeds (Fig. 2). The largest pairwise dif- 
ferences were seen between Spanish Water Dog and Wheaten Terrier (OR = 7.20, p = 0.001), Spanish Water 
Dog and Finnish Lapponian Dog (OR = 6.53, p = 0.001), Spanish Water Dog and Labrador Retriever (OR = 
6.52, p = 0.001), Spanish Water Dog and Jack Russell Terrier (OR = 6.43, p = 0.001), and Shetland Sheepdog 
and Wheaten Terrier (OR = 6.40, p = 0.001). Significant pairwise breed differences are summarised in the 
Supplementary Table S5 and all pairwise breed differences are presented in the Supplementary Dataset. Based on 
previous publications, we made an a priori hypothesis that Chihuahua, Jack Russell Terrier, Lagotto Romagnolo, 
and Shetland Sheepdog would be more fearful than German Shepherd Dog, Golden Retriever, and Labrador 
Retriever. We detected a significant difference between the groups, with the first breed group being more likely 
afraid of strangers than the latter one (OR = 1.65, p = 0.031) (Supplementary Table S4). 
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Figure 2. The effects of breed, sex, sterilisation, socialisation, living environment, and body size on fear of 
strangers in the logistic regression analysis. (a) Breeds differed in their likelihood of showing fear of strangers, 
with Spanish Water Dog being the most and Wheaten Terrier the least fearful breed. (b) Female dogs were more 
afraid of strangers than male dogs. (c) Intact dogs showed less fear of strangers than neutered dogs. (d) Dogs 
that had less socialisation experiences in the age of 7-16 weeks showed more fear of strangers. (e) Dogs living in 
a more urban environment had a higher likelihood of showing fear of strangers. (f) Small sized dogs were more 
afraid of strangers than medium sized dogs. Grey lines (d,e) and error bars (a,b,c,f) indicate 95% confidence 
limits. N = 5,932. 


Dogs showing fear of strangers were more often females (OR = 0.85, p = 0.027) and neutered (OR = 0.62, p < 
0.001) (Table 1, Supplementary Table S4, Fig. 2). Body size was also associated with fear of strangers, but the asso- 
ciation was evident only between small and medium sized dogs (OR = 1.50, p = 0.01) (Supplementary Table S4, 
Fig. 2). Additionally, fearful dogs lived in more urban environments (x? = 14.54, DF = 1, p = 0.001) (Table 1, 
Fig. 2) and participated less frequently in activities and training. Dogs participating in activities only seldom or 
never were more likely fearful than dogs participating sometimes (OR = 1.36, p = 0.013) or weekly (OR = 1.44, 
p = 0.001) (Supplementary Fig. S2). 

Age was not associated with fear of strangers (linear effect: x? = 0.24, DF = 1, p = 0.626; quadratic effect: ,* 
= 3.27, DF = 1, p = 0.070), and the overall effects of family size and weaning age of the dog were not significant 
after controlling for false discovery rate (Table 1, Supplementary Fig. $2). However, when examining pairwise 
comparisons between dogs living in different sized families, we found that dogs living with one (‘single’) or two 
adults (‘couple’) were less likely to show fear of strangers that dogs living in families with two children (single vs. 
two children family: OR = 0.67, p = 0.015, couple vs. two children family: OR = 0.74, p = 0.046) (Supplementary 
Table S4, Supplementary Fig. S2). In addition, dogs weaned later than eight weeks of age had a higher likelihood 
to be afraid of strangers than dogs weaned at a normal age (from seven to eight weeks of age; OR = 0.80, p = 
0.014) (Supplementary Table $4, Supplementary Fig. S2). 


Discussion 

We investigated the demographic and environmental factors associated with two subtraits of canine social fearful- 
ness, fear of dogs and fear of strangers, in a large behavioural survey data of 6,000 dogs using logistic regression. 
Social fearfulness was associated with several factors, which were almost identical in both subtraits. To sum up, 
fearful dogs were less socialised during puppyhood, small in body size, females and neutered, and participated 
less frequently in training and other activities. We also found a novel association between the living environment 
of the dog and social fearfulness, as dogs living in a more urban environment were more likely afraid of dogs and 
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strangers. Moreover, we identified significant behavioural differences between dog breeds. These results reveal a 
comprehensive set of demographic and environmental factors associated with social fearfulness, but also high- 
light the role of genetic predisposition. 

The level of socialisation during puppyhood (between seven weeks and four months of age) had the strongest 
association with both fear of dogs and fear of strangers. Fearful dogs had experienced less socialisation events 
in the age of 7-16 weeks than dogs showing no fear. Many mammals, including dogs, have a sensitive period for 
socialisation in the early postnatal life, during which the nervous system is immature and receptive for novel 
external stimuli**’’. Experiences and events taking place especially during 3-14 weeks of age can significantly 
affect dog’s behaviour throughout life??“°°°. In Finland, puppies are usually weaned and adopted to their new 
homes at the age of 7-8 weeks. This life change takes place during the sensitive period for socialisation, making 
the first weeks in the puppy’s new home crucial for the socialisation process. 

Dogs weaned later than eight weeks of age were more likely to show fear of strangers than dogs weaned at a 
normal 7-8 weeks of age. This result might indicate that dogs weaned later encounter less socialisation experi- 
ences during the sensitive period than puppies moving to their new homes earlier. The breeder might not be able 
to provide an adequate amount of different kind of socialisation experiences when compared to the puppy’s new 
family. However, the association of inadequate socialisation and late weaning age might also reflect fearfulness 
already present at puppyhood. A fearful puppy might have to wait longer before it is adopted, and if it shows fear 
in the new home, it may not often meet unfamiliar people. Nevertheless, our results are consistent with earlier 
research, as it has been shown that early socialisation is necessary for normal social behaviour”””?”> and that 
puppies weaned before the age of 12 weeks are bolder than puppies weaned later*’. Therefore, development of 
social fearfulness may strongly depend on the diversity, amount, and proper timing of socialisation experiences. 

We identified breed differences in fear of dogs and strangers. Spanish Water Dog and Shetland Sheepdog had 
a high probability of fear of dogs and strangers, whereas Wheaten Terrier had a low probability of both subtraits. 
Interestingly, the odds of showing fear greatly differed between the most and least fearful breeds, suggesting that 
some breeds might be more susceptible to fear than others. Many breed differences were similar in both fear of 
dogs and strangers, which could indicate comorbidity between these social fearfulness subtraits. However, some 
breeds only showed fear of dogs or strangers but not both. For example, Chihuahua and Border Collie rated high 
in fear of dogs but not in fear of strangers, whereas Miniature Schnauzer rated high in fear of strangers but not in 
fear of dogs. 

Previous research has identified certain breeds that are more fearful than others**°*. These more fearful 
breeds include Chihuahua, Jack Russell Terrier, Lagotto Romagnolo, and Shetland Sheepdog. In contrast, German 
Shepherd Dog, Golden Retriever, Labrador Retriever, and Staff. Bull Terrier have been identified as less fearful 
breeds. We contrasted these fearful and non-fearful breeds and indeed discovered a difference in the probability 
of fear of dogs and strangers. However, more precise comparisons with previous research are difficult as the breed 
composition of study cohorts has been highly variable and moreover, instead of individual breeds, breed groups 
are often used to overcome small sample sizes in the statistical analyses. As a result, comparisons between findings 
of different studies might not be very meaningful. 

We discovered that small dogs were more likely to show fear of dogs and strangers than large dogs. This 
result agrees with previous studies showing that fear of dogs and strangers decreases as the height of the dogs 
increases*’, and that larger dogs are bolder in general***°“'. Small sized dogs may feel more threatened by people 
and larger dogs due to the larger relative size difference and hence be more fearful. In addition, small dogs may 
not be trained and socialised enough”, as their small size makes them easier to handle even when they behave 
badly or unpredictably. Thus, undesirable behaviours may be more tolerated in small dogs than in larger conspe- 
cifics. For the same reason, behaviour may not be an important factor when making breeding plans, and thus, 
unwanted behaviours may accumulate. Another explanation is that body size could be genetically associated with 
reactivity and stress tolerance in dogs. However, as the most fearful breeds also include small sized breeds, despite 
having size as an explanatory variable in the models, the breed itself may be more important factor than the body 
size. This result emphasises the role of genetic predisposition in the development of behaviour and behavioural 
problems. 

Interestingly, the living environment of the dog was associated with social fearfulness, as dogs living in an 
urban environment were more likely afraid of dogs and strangers than dogs living in more natural and agricul- 
tural areas. The relationship between behaviour and the living environment has not been studied in dogs before, 
but studies in humans have reported higher rates of mental disorders in urban areas when compared to rural 
environments***’. However, the results have been somewhat inconsistent depending on, for example, the geo- 
graphical location studied. In addition, the causality and mechanisms behind these associations are still mainly 
uncovered but studies have suggested several factors, such as diet, environmental toxins, stress, and social iso- 
lation****. As dogs share the environment with us, similar environmental factors could mediate this association 
in both species. In addition, urban areas, such as cities, may be very hectic and stressful environments to live as 
they are full of different stimuli, such as sudden and loud noises, potentially predisposing dogs to fear-related 
and other behavioural problems. On the other hand, dogs living in rural environments may meet other dogs and 
strangers less frequently due to the lower population density in rural areas. Nevertheless, in the lack of previous 
literature, no further conclusions can be made about the observed association between the urban environment 
and higher social fearfulness in dogs. 

We found that dogs engaging in activities and training less frequently were more likely to show fear. In both 
subtraits, dogs participating only seldom in activities and training had a significantly higher likelihood to be 
fearful when compared to dogs participating in activities and training sometimes or weekly. Additionally, fear of 
conspecifics was associated with daily exercise. Dogs exercising less than two hours a day were more likely fearful 
than dogs getting more than three hours of exercise. Interestingly, lower levels of exercise have been previously 
associated with noise sensitivity and fear of strangers””*”. The observed associations may be explained by multiple 
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different phenomena. Activities and exercise can fulfil dogs’ species-specific needs, acting as a way to reduce 
stress*’. Therefore, dogs participating in activities more often may cope better in challenging situations. Previous 
research has demonstrated that both social and environmental enrichment together with exercise improves the 
wellbeing of dogs in kennels*”°°. Moreover, dogs that engage more in activities and exercise may also habituate 
to the fearful stimuli. Additionally, dogs that are trained often likely interact more with their owners, possibly 
strengthening the human-dog bond”*". Thus, the dog may feel safer in the presence of its owner and show less 
fear. On the other hand, owners may not be willing to train and exercise fearful dogs, as they may be difficult to 
handle and behave unpredictably, causing inconvenience and stress to the owner. Therefore, more research is 
needed to better understand the relationship between activities and social fearfulness. 

Consistent sex differences in fearfulness were observed both in our dataset and in previous studies 
Females tend to be generally more fearful than males, and we observed the same trend in both fear of dogs and 
fear of strangers. In addition, neutered dogs were more likely afraid of dogs and strangers than intact dogs. Similar 
results have also been obtained earlier’***°*. Moreover, we detected a significant interaction between sex and 
sterilisation in fear of dogs, as intact males were less fearful than intact females, but neutered males and females 
did not differ from each other. Interestingly, no interaction between sex and sterilisation was observed in fear 
of strangers. In some previous studies, neutering has been demonstrated as a risk factor for several problematic 
behaviours, potentially because of the hormonal changes caused by gonadectomy**”’. However, sometimes vet- 
erinarians may suggest neutering as a treatment for behavioural problems. In addition, dogs with undesirable 
behaviours are less likely used for breeding purposes, and thus they may be more likely neutered. Therefore, the 
causality of the relationship between neutering and social fearfulness remains elusive in our study. 

Older dogs showed less fear of dogs, as the likelihood of fearfulness decreased after six years of age. A similar 
trend was also observed in fear of strangers, but the association was not significant. These results are supported 
by our previous studies: Tiira and Lohi (2015)”* discovered that younger dogs were generally more fearful, and 
Salonen et al.'” found that dogs 4-8 years old were most fearful of both dogs and strangers. However, opposite 
results have also been demonstrated**”’. The decline in fearfulness with age suggests that dogs may habituate to 
fearful stimuli and develop strategies to cope in these situations. Moreover, owners of fearful dogs may also adapt 
their actions, and learn to avoid situations in which the dog feels uncomfortable and shows fear. 

In this survey study, we demonstrate that socialisation during puppyhood is strongly associated with social 
fearfulness in dogs, consistent with previous research. We have also replicated other findings from previous stud- 
ies, as we show that fearful dogs are more often small, females and neutered, and participate less often in training 
and exercising. In addition, we identified several breed differences, suggesting that some breeds may be more vul- 
nerable to develop social fear-related problems than others. Moreover, we report a novel association between the 
living environment of the dog and social fearfulness that requires further research. Most of the identified risk fac- 
tors were common for both fear of dogs and fear of strangers, but we also identified some subtrait-specific risks, 
such as the association of daily exercise and fear of dogs, and the association of weaning age and fear of strangers. 
Our results indicate that social fearfulness in dogs is affected by multiple demographic and environmental factors 
and suggest that careful consideration and management of these factors could improve the wellbeing of pet dogs. 


22,33,36,37,52,53 


Methods 

Data collection. An online owner-filled questionnaire was designed to collect information on canine 
behaviour in a pet dog population’’. The questionnaire was divided into seven main sections based on different 
behavioural traits, including fear of dogs, strangers, and novel situations (labelled as ‘fear’); aggression toward 
human family members and strangers (labelled as ‘aggressiom); fear of thunder, fireworks, and gunshot (labelled 
as ‘noise sensitivity’); fear of surfaces and heights; inattention and hyperactivity/impulsivity; separation anxiety; 
and compulsive behaviour. In addition, the questionnaire included several questions concerning the background 
and the current living environment of the dog, such as socialisation experiences during the age of 7-16 weeks, 
activities that the dog possibly participates in (e.g. dog shows, agility, or herding), amount of daily exercise, and 
the number of other dogs in the family. All the questions included in the behaviour questionnaire can be found as 
Supplementary material in the paper of Salonen et al.”’. 

In the questionnaire section concerning fear, the potential fearful reactions towards strangers and other dogs 
(labelled as ‘social fearfulness’) as well as novel situations (labelled as ‘non-social fearfulness’) were asked. As this 
paper focuses on studying the demographic and environmental factors associated with social fear in dogs, only 
the structure of the social fear subdivision is described in more detail here. If the owner answered that his/her 
dog reacts fearfully when meeting a stranger and/or a strange dog, we required him/her to report how often these 
reactions occur (rarely, 0-20% of occasions; sometimes, 20-40% of occasions; often, 40-60% of occasions; almost 
always, 60-100% of occasions; always, 100% of occasions). In addition, we asked him/her to more specifically 
indicate how the dog behaved in these situations (e.g. the dog withdraws when meeting a stranger). Moreover, if 
the owner responded that his/her dog does not react fearfully in those situations, he/she was asked to describe the 
dog’s behaviour more specifically to give us possibility to evaluate the behavioural reactions of the dog. The fear 
section of the questionnaire was previously shown to have high reliability and validity”. 

As our aim was to reveal demographic and environmental factors associated with both fear of dogs and fear 
of strangers, two separate binomial (event/non-event) variables (‘fear of dogs’ and ‘fear of strangers’) describ- 
ing these behaviours were derived from the questionnaire data. According to the owners, non-fearful dogs are 
reported to indicate no fear of other dogs (‘fear of dogs’) or strangers (‘fear of strangers’) whereas fearful dogs 
are reported to show fear of strange dogs (‘fear of dogs’) or strangers (‘fear of strangers’) in more than 40% of the 
occasions (Supplementary Table S6). 

The questionnaire was advertised to dog owners on Facebook and via breed organisations. We obtained an 
informed consent from all participants and informed them that the questionnaire answers would be used for 
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research. We emphasised that all data will remain confidential and that individual dogs and owners cannot be 
identified from the published results. Questionnaire replies were collected during 2015-2018. 


Demographic and environmental variables. Before statistical analyses, some demographic and envi- 
ronmental variables derived from the questionnaire were edited and some new variables were created. First, only 
breeds with enough individuals in both non-fearful and fearful groups (>10 individuals/group) were selected for 
the analyses. Individuals in other breeds were grouped together under ‘other’ breed group. The selected breeds 
included Border Collie, Cairn Terrier, Chihuahua, Chinese Crested Dog, Coton de Tuléar, Finnish Lapponian 
Dog, German Shepherd Dog, Golden Retriever, Irish Soft Coated Wheaten Terrier (labelled as “Wheaten Terrier’), 
Jack Russell Terrier, Labrador Retriever, Lagotto Romagnolo, Lapponian Herder, Medium size Spitz, Miniature 
Poodle (including Toy, Miniature, and Medium Poodle), Miniature Schnauzer, Pembroke Welsh Corgi, Rough 
Collie, Shetland Sheepdog, Smooth Collie, Spanish Water Dog, and Staffordshire Bull Terrier (labelled as ‘Staff. 
Bull Terrier’) (Supplementary Table S1). Mixed breed dogs were also included in the data. In the ‘fear of strangers’ 
analysis, Cairn Terrier, Pembroke Welsh Corgi, and Staff. Bull Terrier were included in the group ‘other’ as they 
had too few individuals for the logistic regression analysis. 

Second, we created a new categorical variable, body size, based on the average heights of breeds. The heights 
were drawn from FCI and AKC standards, when available. For other breeds, with no FCI or AKC standards, the 
average heights were determined based on heights reported by breed clubs. If different heights were reported 
for females and males, the mean was calculated and used as the average height. Individuals with breed average 
heights of <35.cm, 36-49 cm, or >50cm belonged to small, medium, or large categories, respectively. When body 
size was included in the analyses, mixed breed dogs were excluded as their body size could not be determined. 

Third, a continuous variable labelled as socialisation score was calculated based on the frequency of sociali- 
sation events when the dog was 7-16 weeks old. The score is a sum of the frequencies (0 = never; 1 = 1-2 times 
during the puppyhood; 2 = 1-2 times during the puppyhood to 2 times per month; 3 = twice a month to twice a 
week; 4 = twice a week to once a day; 5 = several times a day) the dog met unfamiliar men, women and children, 
unfamiliar adult dogs, visited city or other place with traffic and many people, and travelled by car or by bus. 

Fourth, a continuous variable labelled as urban environment score describing the environmental land-use 
around the dog’s home was created. We derived the geographical coordinates for each home from addresses pro- 
vided by the dog owners, and calculated the coverage of three land-use types, artificial surfaces, agricultural areas, 
and forests and semi-natural areas within a three-kilometre range around the homes using the land-use database 
CORINE2012. Finally, we simplified the coverages into one continuous variable with principal component anal- 
ysis (PCA). The higher the urban environment score, the more urban was the environment. 

Variables derived from the questionnaire data and included in the analyses are described in more detail in the 
Supplementary Table S6. 


Statistical analyses. Logistic regression was used to find demographic and environmental factors asso- 
ciated with two subtraits of social fearfulness, fear of dogs and fear of strangers. Initially, the data consisted of 
13,715 dogs in 264 breeds. After including only fearful and non-fearful dogs and excluding individuals with 
missing or incomplete responses, the data consisted of 5,343 dogs in ‘fear of dogs’ and 5,858 in ‘fear of strangers. 
The subtraits were used as binary response variables in the analyses: fearful and non-fearful dogs constituted 
the event and the non-event, respectively. Based on previous literature, we included several explanatory vari- 
ables in the analyses. Demographic explanatory variables included age, sex, sterilisation, breed, and body size. 
Environmental explanatory variables included socialisation score, weaning age, urban environment score, daily 
exercise, activities/training, owner’s dog experience, dogs in the family, family size, and daily time spent alone. 
To find the models with the best fit, a forward stepwise model selection by Akaike Information Criterion (AIC) 
values was used and was initiated with a starter model of sex and age in both ‘fear of dogs’ and ‘fear of strangers’ 
analyses. The AIC model selection process and the final models are presented in the Supplementary Table $7. In 
the ‘fear of dogs’ model, the interaction between sex and sterilisation was significant and improved the AIC value 
of the model. Therefore, the interaction term (sex*sterilisation) was included in the final model. To maximise the 
sample sizes, new subsets of the initial data were created to include only those variables included in the final mod- 
els. This resulted in datasets of 5,973 and 5,932 individuals for ‘fear of dogs’ and ‘fear of strangers, respectively. 

After model selection, generalised additive models were fitted with the package ‘gam’ in R to test the linearity 
assumption of continuous variables. If the linearity assumption of a variable was not met, the variable squared was 
also included in the model (e.g. age*2). Packages ‘broom’ and ‘dplyr’” in R were used to inspect possible outliers 
in the datasets. Standardised residuals were plotted using the package ‘ggplot2’® in R. Generalised variance infla- 
tion factor (gVIF) was utilised to test multicollinearity with the package ‘car’*’, and the area under the receiver 
operator characteristic curve (AUC) was calculated to estimate how well the model predicted the event and the 
non-event using the package ‘pROC’™ in R. The predictive abilities of both models were reasonable (AUC(fear of 
dogs) = 0.728 and AUC(fear of strangers) = 0.688). 

Estimated marginal means were calculated for categorical explanatory variables with the package ‘emmeans’®? 
in R. The effects of continuous explanatory variables (adjusting for other variables in the models) were obtained 
with the package ‘effects’ in R. To obtain the overall effects of the explanatory variables, analysis of variance 
(ANOVA) was conducted with the package ‘car’*! in R. Based on previous literature, we determined several con- 
trasts between levels of explanatory variables a priori. First, we requested a contrast between potentially more 
fearful (Chihuahua, Jack Russell Terrier, Lagotto Romagnolo, and Shetland Sheepdog) and less fearful (German 
Shepherd Dog, Golden Retriever, Labrador Retriever, and Staff. Bull Terrier) breeds****. Second, we hypothesised 
that female dogs would be more fearful than male dogs*’**?”*°?, Third, we hypothesised that large dogs would 
differ from small dogs in their behaviour**“'. Fourth, contrast between early weaning (weaned <7 weeks of age) 
and normal weaning age (7-8 weeks of age) as well as contrast between normal weaning age and late weaning age 
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(>8 weeks of age) were requested”*”?. Fifth, we hypothesised that younger dogs would be more fearful”. Finally, 
we made a hypothesis that dogs with less socialisation (i.e. lower socialisation scores) would show more fearful 
behaviour”. The hypotheses and contrasts were the same in both ‘fear of dogs’ and ‘fear of strangers’ analyses, 
except for weaning age which was only included in the ‘fear of strangers’ model. Furthermore, in the ‘fear of 
strangers analysis, the breed contrast was as follows: Chihuahua, Jack Russell Terrier, Lagotto Romagnolo, and 
Shetland Sheepdog versus German Shepherd Dog, Golden Retriever, and Labrador Retriever (Staff. Bull Terrier 
was left out due to a low sample size). 

A priori contrasts and all pairwise comparisons between levels of the included categorical variables were 
examined with the package ‘emmeans’® in R. As we had several categorical variables, the number of pairwise 
comparisons was high and therefore all obtained p-values, except contrasts chosen a priori, were controlled for 
false discovery rate (FRD). The significance cut-off was set at p-value <0.05. 

All statistical analyses were conducted in R version 3.6.1. 


Data availability 


The anonymised data used in this study can be found as Supplementary material in the paper of Salonen et al.!”. 
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